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ABSTRACT 

We identify a population of morphologically defined E/SO galaxies lying on the locus of late-type 
galaxies in color-stellar mass space — the "blue sequence" — at the present epoch. Using three 
samples (from the Nearby Field Galaxy Survey or NFGS, a merged HyperLeda/SDSS/2MASS catalog, 
and the NYU Value Added Galaxy Catalog), we analyze blue-sequence E/SOs with stellar masses >10 8 
Mq, arguing that individual objects may be evolving either up toward the red sequence or down into 
the blue sequence. Blue-sequence E/SO galaxies become more common with decreasing stellar mass, 
comprising <2% of E/SOs near the "shutdown mass" M s ~ 1-2 x 10 11 Mq, increasing to >5% near 
the "bimodality mass" M b ~ 3 x 10 10 Mq, and sharply rising to >20-30% below the "threshold 
mass" M t ~ 4-6 x 10 9 Mq, down to our completeness analysis limit at ^10 9 Mq. The strong 
emergence of blue-sequence E/SOs below M t coincides with a previously reported global increase in 
mean atomic gas fractions below M t for galaxies of all types on both sequences, suggesting that 
the availability of cold gas may be basic to blue-sequence E/SOs' existence. Environmental analysis 
reveals that many sub-Mj, blue-sequence E/SOs reside in low to intermediate density environments. 
Thus the bulk of the population we analyze appears distinct from the generally lower-mass cluster 
dE population; SO morphologies with a range of bulge sizes are typical. In mass-radius and mass- 
a scaling relations, blue-sequence E/SOs are more similar to red-sequence E/SOs than to late- type 
galaxies, but they represent a transitional class. While some of them, especially in the high-mass 
range from Mj, to M s , resemble major- merger remnants that will likely fade onto the red sequence, 
most blue-sequence E/SOs below Mb show signs of disk and/or pseudobulge building, which may be 
enhanced by companion interactions. The blue overall colors of blue-sequence E/SOs are most clearly 
linked to blue outer disks, but also reflect blue centers and more frequent blue-centered color gradients 
than seen in red-sequence E/SOs. Notably, all E/SOs in the NFGS with polar or counterrotating gas lie 
on or near the blue sequence, and most of these systems show signs of secondary stellar disks forming 
in the decoupled gas. From star formation rates and gas fractions, we infer significant recent and 
ongoing morphological transformation in the blue-sequence E/SO population, especially below Mf,. 
We argue that sub- Mb blue-sequence E/SOs occupy a "sweet spot" in stellar mass and concentration, 
with both abundant gas and optimally efficient star formation, which may enable the formation of 
large spiral disks. Our results provide evidence for the importance of disk rebuilding after mergers, 
as predicted by hierarchical models of galaxy formation. 
Subject headings: galaxies: evolution 



1. INTRODUCTION 

Modern galaxy surveys find that most galaxies occupy 
two distinct loci in color-stellar mass space, the "red se- 
quence" and th e "blue sequence" ( or "cloud" or "dis- 
tribution," etc.; iStrateva et al J 120011 iBaldrv et all 12004 
iBell et al]l2004j ). Because the high-mass end of the red 
sequence corresponds to the well-known color-magnitude 
relation of cluster E/SO galaxies, it is natural to identify 
red-sequence galaxies with E/SO types and blue-sequence 
galaxies with spiral/irregular types, except for some con- 
tamination of the red sequence by dust-reddened late- 
type systems. 

Here we show that this basic morphology-color corre- 
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spondence begins to fail for stellar masses M* < 1-2 x 
10 11 Mq, as morphologically defined E/SO galaxies start 
to appear on the blu e sequence in z = surveys (see also 
iBamford et al.ll2009ft . The failure becomes dramatic be- 
low a threshold mass M t of 4-6 x 10 9 Mq, matching the 
mass scale below which the mean atomic gas content of 
low-z galaxies increases substantially on bo th sequences 
()Kannappanll2004l : iKannappan fc Weill2008t ). 4 

A specific stellar mass scale linked t o shifts in galaxy 
prope rties was first highlighted by Kau ffmann et al.l 
(2003), who found that galaxy star formation histo- 
ries (SFHs) are qualitatively different on either side of 
a crossover or bimodality mass Mb ~ 3 x 10 10 Mq. 
Their analysis shows the transitional mass range start- 

4 Kannappan (2004) incorrectly associated the gas-richness 
th reshold mass with the 3 X 10 10 Mq "bimodality" mass scale 
of IKauffmann e~ al. (2003) based on an assumed correspondence 
between t h e stel lar mass calibrations of Kauffmann et al. and 
IBell et aLl (120031). On f urther investigation (see both § 12.31 and 
Kannappan & Wei 2008), the gas-richness threshold mass in fact 
corresponds to ~4-6 X 10 9 Mq and marks the lower edge of a 
transitional mass range centered on the bimodality mass. 
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ing at log M*/M ~ 9.5-10, below which galaxies have 
bursty SFHs, and ending at log M*/M ~ 11-11.5, be- 
yond which SFHs reflect uniformly ancient stellar pop- 
ulations. Structural properties also change, with high- 
mass galaxies having higher concentration, higher sur- 
face brightness, and a higher ratio of spher oids to disks 
(|Kauffmann et all 120031 (Driver et al"1l2006f). In a more 
focuse d study of edge-on disk galaxies. iDalcanton et al.1 
( 2004) report the emergence of bulges as well as efficient 
star formation, as evidenced by prominent dust lanes, 
above V c ~ 120 km s _1 (corresponding to log M*/M ~ 
9.5 in t he stellar-mass Tul ly-Fisher relation, see § 13.21 
below) . iBaldrv et alJ (|2004l ) tie mass-dependent shifts in 
galaxy properties to the red and blue sequences by show- 
ing that the red sequence becomes numerically dominant 
above M* - 2-3 x 10 10 M Q . 

Just about every idea advanced to explain why 
galaxy properties change at specific mass scales in- 
volves gas physics. Abrupt shifts in star formation 
efficiency and gas richness near M t have been linked 
to the interplay of gas infall, superno va-driven winds, 
and c hanges in surface mass density (|Dalcanton et al.l 
120041: IDalcanton! l2007h . The buildup of the red se- 
quence above Mb has been linked to "quenching" pro- 
cesses that shut down star formation, such as gas 
loss caused by AGN feedback, consumption of all 
available gas in the wake of violent merging, ram- 
pressure stripping/harassment/strangulation in clusters, 
and/or shock- heating that turns off efficient "cold- 
mode " gas accretion (IBaldrv et al.l 120041: iKennev et all 
2004 ISpringel et all I2005t iDekel fc Birnboiml 120061: 
Cattaneo et al.1 120061 : iFaber et all l2007f ). Some combi- 



nation of these processes may be required to ensure that 
quenching is efficient, p ermanent, and mass-dependent 
(jDekel fc Birnboim|[200l) . 

Higher-redshift studies now suggest that the crossover 
mass between disks and spheroids has evolved down- 
ward over time from M* ~ 1-2 x 10 11 M at z ~ 1 
- notably, equal to the shutdown mass M s marking 
the upper end of the transitional m ass range today - 
to today's crossover ma ss near Mh_(Bundv e t al.l 120051 : 
iFranceschini et all 120061 : see also ICimatti et al.l |2006| ). 
This tre nd is coincident wi th many other signs of "down- 
sizing" (jCowie et al.lll996l ) in the galaxy formation pro- 
cess, including the downwardly evolving upper mass 
thres h old for strong star formation (e.g . . lJuneau et al.l 
120051: iKodama et all 12004 iTresse et al.l |2007D as well 
as a similarly evolving upper threshold for irregular 
morphologies, abundant blue co mpact galaxies, and 
strong kinematic disturbances (|Mallen-Ornelas et al.l 
ll999tlKannappan fc Bartonll2004h . Whether downsizing 
can be traced back to still higher redshifts is unclear, 
as the red sequenc e is only tentative ly detected beyond 
z ~ 1 .5 (comparing lKriek et al.ll200"8l and lCirasuolo et al.l 
120071 ). Semianalytic models that include cold accretion 
predict that the highest mass galaxies, residing in ~ 
10 12 Mq dark matter halos (implying stellar masses near 
M s ), formed rapidly at very high redshift before devel- 
oping stable shocks th at inhibited further star formation 
(|Cattaneo et al.ll2006f h 

From these results, we infer the importance of three 
mass regimes for our z — analysis: (1) near the 
shutdown mass M s , at which blue-sequence E/SOs first 
emerge, and above which nearly all galaxies are old and 



red; (2) between M s and M t , where the signatures of 
downsizing from z = 1 — may remain; and (3) below the 
threshold mass M t , where blue-sequence E/SOs, bursty 
SFHs, and high gas fractions become suddenly abundant. 
The surveys considered here include few galaxies below 
M* ~ a few x 10 8 M , so we largely omit the abun- 
dant cluster dE population that dominates the faint end 
of the red sequence. In fact, many E/SOs between 10 s - 
lO 10 5 Mq o ccupy low-to-intermediate density envir on- 
ments (§ I3TT1 iHogg et al.ll2003t [Trentham et~atll2005l ). 

Given the close association of blue-sequence E/SOs 
with mass regimes characterized by active or recent 
evolution, their intriguing mismatched color-morphology 
status could indicate a transitional state in color (e.g., 
a fading starburst galaxy after a merger), morphology 
(e.g., a disk-building system), or both (e.g, a gas-rich 
merger remnant regenerating a young disk from tidal de- 
bris). Put differently, blue-sequence E/SOs may be evolv- 
ing from blue to red, from red to blue, or wholly within 
the blue sequence. Or they may not be evolving much at 
all. 

In this work we document the blue-sequence E/S0 pop- 
ulation systematically for the first time, exploring prop- 
erties and demographics to constrain the disk-building 
and fading-starburst scenarios. We ask the reader to 
suspend all preconceptions about E/S0 colors, gas con- 
tents, masses, environments, and even bulge-to-disk ra- 
tios, as many accepted "facts" about E/SOs reflect the 
properties of the high-mass red sequence. Section [2] de- 
scribes our survey samples, our morphological definition 
of an E/S0, how we measure stellar mass, and how we 
divide the red and blue sequences. Section [3] provides a 
demographic overview of the blue-sequence E/S0 popula- 
tion, including frequency among all E/SOs, distribution 
by mass and environment, and dynamical scaling rela- 
tions. Section |4] examines the detailed characteristics of 
blue-sequence E/SOs, including morphological substruc- 
ture, gas content and star formation, and clues to in- 
teraction status. Section [5] follows up on an intriguing 
association of gas-stellar counterrotation and polar rings 
with blue-sequence E/SOs, providing an in-depth analy- 
sis of new stellar kinematic data for five such galaxies to 
probe the possibility of secondary stellar disk growth in 
these systems. Finally, Section [6] reviews our results in 
the context of different evolutionary scenarios and dis- 
cusses implications for disk formation from z = 1 — 0. 

2. DATA AND METHODS 

We assume Hq = 70 km s _1 Mpc _1 and d — cz/Hq 
throughout. At the low redshifts of our sample galaxies, 
neglecting A introduces errors <0.02 dex in stellar mass. 

2.1. Samples 

We analyze two primary samples, one drawn from 
the Ne arby Field Galaxy Survey (NFGS, Uansen et al.l 
l2000bl hereafter J00) and the other drawn from the 
HyperLeda, Sloan Digital Sky Survey (SDSS) DR4, 
and Two Micron All-Sky Survey ( 2MASS) databases 



(hereafter, "HyperLeda+" samp le: iPaturel et~a 



lAdelman-M cCarth v" et all 120061 : Uarrett et al 



2003a; 
2000). 



These samples are shown in Fig. [TJ The NFGS sam- 
ple offers a homogeneous and high-quality data set while 
the HyperLeda+ sample offers the statistical advan- 
tage of thousands of galaxies. A third sample, the 
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Fig. 1. — Distribution of morphologies in color vs. stellar mass space, for (a) the HyperLeda+ sample and (b) the NFGS sam ple. In both 
panels, the dashed line divides the red sequence from the blue sequence, wit h a c alibrated offset from u — r to U — R (§ 12.41 1 and a shift 
of 0.04 dex between m asses determined with SDSS vs. NFGS photometry (§ 12. 31 1. The gray curve in the panel a is the sequence divider 
of IBaldrv et al.l ([2004), approximately converted from L r to M„ units. Blue-sequence E/SOs in the HyperLeda+ sample have been vetted 
as described in § 12,11 but late types and red-sequence E/SOs in the HyperLeda-l- sample have not been reclassified; thus the blue-sequence 
E/SOs shown in panel a represent a lower limit on the fractional size of the population. Lighter gray symbols in panel b show galaxies with 
visually noted defects in 2MASS photometry; see footnote [6] 



NYU Val ue-Added Galaxy C atalog low-redshift sample 
(VAGC, iBlanton et al.l l2005h . is used to evaluate com- 
pleteness in the HyperLeda-l- sample. We describe the 
samples here and amplify on the morphological criteria 
used to identify E/SOs in S [2T2l 

2.1.1. The Nearby Field Galaxy Survey (NFGS) Sample 

The NFGS database includes ^200 galaxies selected to 
fairly represent the natural abundance of morphologies 
in the local universe over a wide range in luminosity. 
Notwithstanding its name, the survey spans a variety 
of environments including clusters. NFGS galaxies were 
given numerical morphological type s as part of the CfA 1 
redshift survey (jHuchra et al.lll983h and later reclassified 
by J00 on the same deVaucoleurs-based system, using 
new high-quality CCD imaging. 5 Because the NFGS was 
selected without explicit bias in color or morphology, the 
properties of blue-sequence E/SOs in this survey should 
be reasonably representative, within the constraints of 
_B-band selection and small number statistics. 

The NFGS sample allows us to examine a broad sample 
of E/SOs with supporting data not available for the Hy- 
perLeda + sample, most notab ly, integrated spectropho- 
tometry ^Jjansen^^i^|20(^^ 

matics (IKannappan k. Fabric ant 200ll lKannappan et al.1 
I21MI2009T see also new data in § 15. 2[) . Other data used 
in analyzing the NFGS sample include UBR photome- 
try and Virgocentric-infall correcte d redshifts from J00 , 
2MASS AT-band photometry from I Jarrett et all (|2000f ). 
and atomic g as masses deri v ed from the homogenized 
HI catalog of iPaturel et all (|2003bf) using Affji+Hc — 
1.4(2.36 x 10 5 /d 2 ) M , where / is the integrated line 
flux in Jy km s , d is the distance in Mpc, and the fac- 
tor of 1.4 represents the helium mass correction. SDSS 
ugriz data are available for ~60% of the NFGS; we use 
these magnitudes to calibrate shifts in colors and mass 

5 We convert the type Pec galaxy UGC9562, a polar ring 
galaxy with an SO host, to type SO for this paper. We also re- 
classify UGC6570 (=A 11332+3536) from SO to SO/a following 
IKannap pan et al. (2004|). In all other cases we use the types given 
by J00 for the NFGS. 



estimates between the U BRJHK (primary NFGS) and 
ugrizJHK (HyperLeda-l-) systems. 

Our primary NFGS sample consists of 164 galaxies, 
including 52 E/SOs. We reject galaxies with missing 
spectra or point-source morphologies. We also exclude 
a possible blue-sequence E/SO, NGC 2824, with severe 
contamination from a nearby bright star; however, we 
retain one red-sequence E/SO, NGC 3605, whose B and 
i?-band magnitudes are mildly affected by blending with 
a companion, but whose U — R color can still be esti- 
mated reasonably well based on stellar population fits to 
the full SED. All sources in this primary sample have 
.ftT-band fluxes with error <0.3 mag. 6 

Both spectra and photometry are corr ected for fore- 
groun d extinction using the maps of ISchlegel et al.l 
(|1998h and, f or the spec t ra, th e Milky Way extinc- 
tion c urve of lO'Donnelj (| 19941 ) as given by iMcCall 
(|200l . We use total, extrapolated magnitudes for stel- 
lar mass estimation (along with spectra, § I2.3|) . The red 
and blue sequences are divided using isophotal U — R 
colors measured within the S-band 26 mag arcsec -2 
isophote, to minimize extrapolation errors. We esti- 
mate fc-corrections from the SED fits used to determine 
stellar masses (see § 12. 3|) . The notation (U — R) 1 sig- 
nifies that we also correct for internal extinction in a 
population-averaged way, by applying the inclination- 
based corrections of iTullv et al.1 (Il998f) adapted to the 
NFG S and SDSS passbands (see IKannappan et al.ll2002l . 
2009). These corrections are appropriate for galaxies 
with gas; for galaxies lacking extended emission lines, 
we assume zero internal extinction. 

6 Detailed inspection reveals galaxies for which 2MASS data are 
available but pipeline data products are defective, due to poor de- 
blendi ng, electronic glitches, c ut-off edges, or other technical prob- 
lems (Kannappan et al. 2009). This flagging would reject some 
galaxies that pass our .ff-band error criterion, but we allow these 
galaxies in this paper to facilitate uniform comparison with the 
HyperLeda+ and NYU VAGC samples. Fig. [Tp indicates galaxies 
that pass the iC-band criterion but not the more detailed inspection 
in a lighter shade of gray; all but one are late- type galaxies. 
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2.1.2. The HyperLeda/SDSS/2MASS (HyperLeda+) and 
NYU V 'AG 'C Samples 

The HyperLeda+ sample is a larger, less uniform data 
set consisting of 3783 galaxies with morphological types 
listed in HyperLeda plus ugr and K magnitudes avail- 
able from the SDSS and 2MASS galaxy catalogs (within 
a 6" matching radius and with magnitude errors < 0.3 
in K and < 0.15 in ugr). We include only galaxies with 
diameter D25 less than 120" to minimize pipeline reduc- 
tion errors in 2MASS and SDSS. To improve classifica- 
tion reliability, we exclude galaxies listed as "multiple" 
in HyperLeda, and we require diameter >40" and ax- 
ial ratio <0.45 (inclination less than ^72°). We adopt 
Virgocentric-infall corrected redshifts from HyperLeda 
when available and otherwise use SDSS redshifts, which 
we correct to the Local Group frame of reference. Data 
for a handful of catalog galaxies that pass our diameter 
cuts but have nominal redshifts cz > 15000 km s" 1 are 
rejected as spurious. We further require cz > 500 km s _1 
and exclude sources with cz < 1500 km s~ x if they lack 
Virgocentric infall corrections. 

For our working HyperLeda+ sample we have re- 
classified all candidate blue-sequence E/S0 galaxies us- 
ing SDSS g-band imaging, with reclassifications homog- 
enized between two independent observers (SJK and 
AJB), taking the NFGS sample as a reference. Starting 
with HyperLeda numerical types <0 (SO/a), our reclassi- 
fication discards ~60% of candidate blue-sequence E/SOs 
as actually having type Sa or later. The remaining con- 
firmed blue-sequence E/SOs represent a lower limit to the 
true population, because we have not retyped galaxies 
initially classified as late type to find mistyped E/SOs. 7 

In analyzing the HyperLeda+ sample we employ to- 
tal magnitudes determined from extrapolated profile fits, 
i.e., SDSS model magnitudes and 2MASS extrapolated 
magnitudes. We determine /c-corrections and foreground 
extinction corrections as for the NFGS sample. Inter- 
nal extinction corrections are also matched in princi- 
ple. However, because we do not know which galax- 
ies have extended ionized gas emission lines (the crite- 
rion used to decide whether to apply an internal ex- 
tinction correction for an NFGS galaxy), we exploit a 
division in mass-radius space that predicts surprisingly 
well which galaxies in the NFGS have extended emission. 
Thus an internal extinction correction is applied when- 
ever the Petrosian r-band 90% light radius 7-90 satisfies 
logrgo/kpc > -4.55 + 0.50 log M*/M©. 

We determine colors and masses for the NYU VAGC 
low-redshift sample almost exactly as for the Hyper- 
Leda+ sample, to facilitate the comparative analysis in 
§ 13.11 However, axial ratios for VAGC galaxies are avail- 
able only by cross-matching with the online HyperLeda 
database, which yields values for only ^60-70% of the 

7 The HyperLeda+ sample used in this paper was defined prior 
to a 2006 update to the HyperLeda database that corrected a sys- 
tematic error in galaxy radii and added new data. Adopting the 
revised radii would slightly change our existing sample, because two 
galaxies would just miss the lower diameter limit, but would not 
otherwise affect our analysis, as the radii we analyze come from the 
2MASS catalog. Using the revised database would also add and re- 
move some candidate E/SOs, because average numerical types can 
shift as new data are added. However, as all of our blue-sequence 
E/SOs are independently verified by inspection of SDSS data, our 
data set remains a robust minimum sampling of the blue-sequence 
E/S0 population. We therefore retain the original sample. 



VAGC galaxies we analyze (depending on sample defini- 
tion, § 13- 1() . Galaxies lacking axial ratio data are treated 
as having no internal extinction. 

2.2. Definition of an E/S0 

Our definition of an E/S0 galaxy is purely mor- 
phological, based on traditional by-eye classification of 
monochrome B or g band images by multiple observers. 
As discussed in § I2.1( we calibrate our reclassifications 
of candidate HyperLeda+ blue-sequence E/SOs against 
the classifications given for the NFGS in J00, which were 
determined on a simplified version of the deVaucoleurs- 
based system used for the CfA 1 redshift survey. Numer- 
ical types < (SO/a) qualify as E/SOs, including all type 
E, cE, SO, and SO/a galaxies. We also include some pe- 
culiar galaxies with predominantly spheroidal structure. 
Our classification system does not identify "blue com- 
pact dwarfs" (BCDs) as a distinct category, but rather 
distributes them among early and late types according 
to their primary morphology. Dwar f ellipticals (as the 
term is used by iBinggeli et alj |1988[ ) fall mostly below 
our mass range, and their tendency to occur in clusters 
is inconsistent with the available environmental data for 
our galaxies (§ 13. 1| . However, we make no explicit ex- 
clusion based on concentration or surface brightness, and 
we find a continuous range of E/S0 properties. 8 

In practice, most E/SOs in the mass range of our anal- 
ysis are SO and SO/a galaxies. Our classification sys- 
tem downplays the importance of bulge-to-disk ratio as a 
defining characteristic of SOs, requiring only a bulge plus 
a smooth outer disk. From this point of view, the main 
criterion distinguishing an SO from a spiral is the pres- 
ence of extended spiral structure in the l atter (as in the 
paral lel sequences classification system of Ivan den Berghl 
1976). Incipient central spiral structure or faint tidal fea- 
tures around a smooth outer disk may lead to an SO/a 
classification. The sharp boundary between SO/a and Sa 
is of course artificial, particularly in the context of our 
interpretation of blue-sequence E/SOs as a transitional 
population. Fig [2] shows three g-band image stretches 
for a galaxy we judge to be literally on the borderline 
between SO/a and Sa. We choose to retain such galaxies 
on the SO/a side of the division; they make up <10% of 
confirmed blue-sequence E/SOs in the HyperLeda+ sam- 
ple. 

Nothing in the morphological definition of an E/S0 
prohibits rings, bars, gas, dust, star formation, or com- 
panion interactions, provided these do not strongly af- 
fect the primary bulge+smooth outer disk morphology. 
Although we largely exclude major mergers in progress 
by rejecting HyperLeda+ systems marked as "multiple," 
our samples include peculiar E/SOs such as polar ring 
systems and spheroids with shells or tails suggesting late- 
stage mergers. 

Fig. [3] shows a representative montage of E/S0 images 
from both sequences and both samples. We stress that 

8 It is worth recalling that the strong dichotomy in dE vs. giant E 
properties found in the Kormendy Relation and other scaling laws 
does not consider SO galaxies, which are typically either excluded 
or included only in their bulge component. SOs display a broad 
range of intermediate morphologies and masses, and SOs formed 
in mergers are likely to exhibit a wide variety of surface bright- 
ness profiles as a function of progenitor gas richness. Therefore an 
inclusive approach is physically motivated. 
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Fig. 2. — Three g-band image stretches for a galaxy on the borderline between SO/a and Sa, UGC4902, which our merged classification 
lists as "X," meaning literally too difficult to label either way. The finger of asymmetric dust to the south and the ambiguous, possibly 
one-sided arm structure make this a challenging galaxy to classify. Such borderline cases make up <10% of confirmed blue-sequence E/SOs 
in the HyperLeda+ sample. 



the color SDSS cutout images in Fig.[3]do not necessarily 
reflect the morphology seen in the monochrome classifi- 
cation images: the cutouts highlight dust and color struc- 
ture while downplaying low surface brightness outer disk 
features. 

Three of the 15 NFGS blue-sequence E/SOs are also in 
the HyperLeda+ sample, but are not considered E/SOs 
in the latter sample. One is UGC9562, the polar ring 
galaxy shown near (logM*, color) = (9.2, 1.3) in Fig. [3J 
because this galaxy is classified as an Sd in the Hyper- 
Leda database, it was never a candidate for reclassifica- 
tion. The other two are UGC6655 and IC1144, which 
are listed as E/SOs in the HyperLeda+ database but 
are among the 60% of candidate blue-sequence E/SOs 
that we rejected upon inspection of SDSS images. Rein- 
spection of the NFGS and SDSS images confirms that 
both are borderline cases, with the NFGS images making 
them look more like SO/a and the SDSS images making 
them look marginally later in type. To ensure uniform 
analysis, we do not change their classifications in either 
sample (i.e., they are plotted as late types in the Hy- 
perLeda+ sample). UGC6655 is shown in Fig. [3] near 
(log A/*, color) = (8.3, 1.5). These discrepancies high- 
light the inexact nature of morphological classification; 
however, § 13.21 will show that E/SOs in both samples 
have structural properties that are quantitatively distinct 
from later type galaxies. 



2.3. Stellar Mass Estimation 

Stellar masses are estimated using an improved ver- 
sion of the stellar population fitting code described in 
iKannappan fe Gaw iser (200^ . The code fits photometry 
(UBRJHK or ugrizJHK) and integrated spectropho- 
tometry (optional: NFGS only) with a grid of models 
b uilt up from the simp le stellar po pulation ( SSP) models 
of iBruzual fe Charlotj (120031) for a lSalpeter! (fl95l Initial 
Mass Function (IMF). We rescale the final masses by a 
factor of 0.7 to approximate an I MF with fewer lo w-mass 
stars (the "diet" Salpeter IMF of lBell et al.1l2003h . Each 
model combines an old SSP (age 1.4, 2.5, 3.5, 4.5... 13.5 
Gyr) plus a young SSP (age 25, 100, 290, or 1000 Myr). 
The young SSP makes up 0%, 1%, 2%, 4%, 8%, 16%, 
32%, or 64% by mass; pure young SSPs are also allowed. 
Each SSP has metallicity Z = 0.008, 0.02, or 0.05 and 
the young SSP can take on 11 extinction/reddening val- 
ue s (ry = 0, 0.12, 0.24... 1-2). F urther details are given 
in IKannappan fe Gawiser] (|2007[) . 

Following iBundv et alJ (|2005h . we define the stellar 
mass estimate not by the best fit, but by the median 
and 68% confidence interval of the mass likelihood dis- 
tribution binned over the grid of models. The distri- 



butions imply typical uncertainties of 0.1-0.2 dex. 9 In 
the binned likelihood approach, the density of models 
can affect results. Our curre nt model grid improves on 
IKannappan fe Gaw iser ( 200?]) in that the old population 
ages are equally spaced through time, and we weight the 
model space such that each of the four young ages con- 
tributes 0.25 x its model likelihood, with the 11 ry's 
further subdividing the likelihood contribution. Note 
that having included ry m the models, we do not ap- 
ply inclination-based internal extinction corrections to 
the input photometry used for mass estimation; how- 
ever, the output likelihood distributions do not strongly 
constrain ry, so inclination-based extinction corrections 
r emain the best option fo r comp uting (U — R) 1 (§ [2]). 

Kannapp an fe Gawiserl (|2007h report that systematic 
errors between different methods of stellar mass estima- 
tion can be as large as factors of 2-3, even with matched 
IMFs. Ou r mass normaliza t ion is fortuitously similar 
to that of iKauffmann et alJ (|2003h . as seen in Fig. |U 
Therefore the bimodality mass identified by Kauffmann 
et al. at 3 x 10 10 M Q has roughly the same value in this 
paper. This absolute mass scale is also similar to that 
found using the original Kannappan & Gawiser code with 
Bruzual-Charlot models, and as such is intermediate be- 
tween t he lower and high er mass scales found using ei- 
ther the Maraston (2005) models or the color-M/L cal- 
ibration oTlBeIOLi2~(2003|), respectively. Thus our ab- 
solute scale is a reasonable compromise. However, the 
relative position of red- and blue-sequence E/SOs vs. spi- 
rals in M*-r and M*-a scaling relations may be sensitive 
to systematics in mass estimation between galaxy classes, 
where relative zero points between classes may vary with 
estimation method, and the possible effects of such sys- 
tematics are discussed in 



2.4. The Red/Blue Sequence Dividing Line 

As shown in Fig. [TJ we divide the red and blue se- 
quences with a line that closely hugs the locus of late-type 
galaxies, leveling at (u — r) 1 ~ 1.8 and 2.3 i n agre ement 
with the functional divider of iBaldrv et al.l ((2004). The 
NFGS line is shifted to the 0.66 mag bluer (U — R) 1 scale 
and the 0.04 dex lower mass scale that result from us- 
ing UBRJHK+spectra, rather than ugrizJHK as input 
data, based on empirical calibration using NFGS galaxies 
with SDSS photometry. We measure the 0.66 mag off- 
set from E/SOs alone to best align the sequence division 
between samples; including Sp/Irr types gives a slightly 
larger 0.75 mag color offset. 

9 Distance errors are added in quadrature after mass estimation, 
as they are the same for all passbands and therefore excluded from 
the fits. 




Fig. 3. — SDSS ugriz color-composite cutouts of E/SOs from both sequences, arranged by u — r color and s tella r mass. We include 
examples from both the HyperLeda+ and NFGS samples, using the u — r to U — R conversion given in Section 12.41 Images are slightly 
shifted to minimize overlap and sized arbitrarily. 
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Fig. 4. — Comparison of s tellar mass estim ates from our 
ugrizJHK fitting code and from Kauffmann ct al. (2003), for 658 
galaxies in common between the HyperLeda-l- and Kauffmann et 
al. samples. The normalization is very similar. 

Additional E/SO systems, possibly analogous to blue- 
sequence E/SOs, lie in the minimally populated zone be- 
tween the sequences (hereafter, "mid-sequence"). We 
separately track the two mid-sequence NFGS E/SOs seen 
in Fig. [TJd near (logM*, color) = (9.4, 1.3), while for the 
HyperLeda+ sample, mid-sequence E/SOs are grouped 
with red sequence E/SOs. 

3. DEMOGRAPHICS AND SCALING RELATIONS 
OF BLUE-SEQUENCE E/SOS 

Having robustly identified a population of blue- 
sequence E/SOs in two independent surveys, we now ex- 
amine how this population compares with red-sequence 
E/SOs and late types in the overall galaxy population. 

3.1. Numbers, Masses, and Environments 

Raw NFGS numbers give the impression that blue- 
sequence E/SOs become dramatically more abun- 
dant with decreasing mass (Fig. [T|), especially below 
logMi,/ Mq ~ 10.5. Here we take a closer look at these 
galaxies' mass and environment distributions. 

We compare the HyperLeda+ sample with the NYU 
VAGC in this analysis, where both are needed because 
the VAGC has a simple selection function without mor- 
phological type data, while the HyperLeda+ sample 
has morphological type data without a simple selection 
function. Fortunately, the HyperLeda+ sample closely 
resembles an apparent magnitude limited (r < 14.9) 
subsample of the VAGC, within similar redshift limits 
(Fig. [5^ and Figs. [5^L-b). 10 Therefore it is reasonable 
to treat the HyperLeda+ sample as an approximate sta- 
tistical sample, and we can correct for color and lumi- 
nosity biases by binning the HyperLeda+ and VAGC 

10 Although the VAGC was designed to improve on the raw 
SDSS, it does not restore large angular-size galaxies shredded by 
the SDSS photometric pipeline and therefore omitted from the 
SDSS redshift survey. This fact presumably explains why Fig. [Bp 
contains more high-mass galaxies than[B^. However, the more dom- 
inant blue sequence in the HyperLeda+ sample likely reflects the 
long history of _B-selected redshift surveys feeding into the Hyper- 
Leda+ database. 



data into mass bins within each color sequence, then 
computing a correction factor for each bin, equaling the 
ratio between the HyperLeda+ counts and the corre- 
sponding VAGC counts for a volume-limited VAGC sub- 
sample (Figs. [5b and [5b; defined by cz = 1000-7500 
km s~ 1 and M r < —17.15, where the latter equals the 
VAGC's r = 18 limit at 7500 km s" 1 ). This volume- 
limited subsample is largely complete in stellar mass to 
a limit of log M*/M Q = 9, as seen in Fig. [51 so we com- 
pute number statistics to this limit. 

Fig. [7] presents both the raw HyperLeda+ red- 
and blue-sequence E/SO counts and the completeness- 
corrected frequency of blue-sequence E/SOs as a function 
of stellar mass. The fraction of E/SOs on the blue se- 
quence rises steadily with decreasing mass, from ^2% at 
M s to ~6% at M b , then shoots up to 20-30% below M t . 
The overall abundance of blue-sequence E/SOs relative 
to the galaxy population rises in tandem, from <0.5% 
to ^2% to ~5%. We have compared results using com- 
pleteness corrections based on the simple volume-limited 
VAGC sample vs. based on a modified volume-limited 
subsample restricted by the same if -band flux and error 
requirements used for the HyperLeda+ and NFGS sam- 
ples. Although the AT-band restrictions reduce overall 
galaxy numbers below log M*/Mq ~ 9.5 (Fig. [Bfc), they 
have minimal effect on relative numbers between the two 
sequences. For consistency, we present percentages com- 
puted with the if-band restrictions. The decline in blue- 
sequence E/SO frequency in the lowest-mass bin of Fig. [7] 
should not be overinterpreted, as shredding and explicit 
surface brightness selection criteria start to significantly 
affect SDSS and thus VAGC dwarf galaxy counts in this 
mass bin. 

We emphasize that all of these blue-sequence E/SO fre- 
quency estimates are lower limits, both because of our 
one-sided morphological reclassification effort (§ 12. 2|) and 
because mid-sequence E/SOs are grouped with the red 
sequence in the HyperLeda+ sample (§ I2.4|) . The un- 
corrected NFGS, which behaves somewhat like the K- 
band restricted version of the VAGC volume- limited sub- 
sample (Fig. EH), shows higher blue-sequence E/SO frac- 
tions. However, these numbers may be affected by the 
survey's i?-band selection, and the small number statis- 
tics and complex selection function of the NFGS make 
completeness-correcting it prohibitive. 

Interestingly, Fig. [7] hints at a connection between red- 
and blue-sequence E/SOs below Mf,: while the two E/SO 
families have very different mass distributions over all 
masses (K-S test probability 2 x 10~ 12 of being drawn 
from the same population), their mass distributions be- 
low Mb are indistinguishable (K-S test probability 0.9 
of being drawn from the same population) . We suggest 
that evolutionary processes affecting both red- and blue- 
sequence E/SOs may be changing across this mass scale, 
plausibly linking the emergence of large numbers of blue- 
sequence E/SOs to the onset of high cold gas fractions, 
bursty star formation, and disky morphologies, starting 
around Mh and becom i ng more pronoun ced below M t 
(|Kauffmann et al.ll2003llKannappanll2004[ ). 

This suggestion agrees with the available environ- 
mental data for red- and blue-sequence E/SOs in the 
logM-t/ Mq = 9-10.5 mass regime, where most E/SOs 
seem to occupy similar, and notably low-density, envi- 
ronments (Fig. [8|) . Completeness corrections as a func- 
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tion of environment are infeasible for either of our sam- 
ples, but within these data sets, the low-density envi- 
ronments of blue-sequence E/SOs appear correlated with 
their intermediate ma sses. This result fi nds support in 
the SDSS analysis of iHogg et al.l ((2003) , who report a 
non-monotonic density trend along the red sequence, 
seen as a dip in typical density for intermediate-mass red- 
sequence galaxies despite higher densities at higher and 
lower masses (the latter likely due to the rise of dEs be- 
low 10 8 -10 9 Mq). Contrary to common assumption, the 
environments of intermediate- mass E/SOs are not neces- 
sarily more dense than those of spira l galaxies, as seen 
in both Fig. 151 and fHogg et al.l <|2003h . Neither result is 
immune to cosmic variance, but both demonstrate that 
intermediate-mass E/SOs are abundant in low-density en- 
vironments where disk regrowth is plausible. Analysis of 
group dynamics would be beneficial for understanding 
the balance of cold gas accretion vs. quenching processes 
in blue-sequence E/SO environments, but such an anal- 
ysis is beyond the scope of this paper. We do consider 
evidence that interactions play a role in blue-sequence 
E/SO evolution in § O and § O] 

Taken together, Figs. [7] and [8] may explain why blue- 
sequence E/SOs have been missed as an important pop- 
ulation up to now, despite their large numbers. First, 
they emerge in numbers competitive with red-sequence 
E/SOs only below M t , so inevitably they are underrep- 
resented in magnitude-limited samples. Second, blue- 
sequence E/SOs are rare in the dense cluster environ- 
ments often targeted for studies of early-type galaxies. 

3.2. M*-r and M„-a Scaling Relations 

Blue-sequence E/SOs are more similar to red-sequence 
E/SOs than to late-type galaxies in both the Af*-radius 
relation and the M*-er (stellar velocity dispersion) rela- 
tion (Figs. [91 and ITU]) . This basic result validates their 
morphological classification. At the same time, both re- 
lations show hints that blue-sequence E/SOs may rep- 
resent a transitional class, and that E/SOs on both se- 
quences change in structure below Mb~M t . The two 
NFGS E/SOs from the mid-sequence zone between the 
red and blue sequences appear most consistent with blue- 
sequence E/SOs in the M*-r and M*-tr relations, but 
these mid-sequence systems fall within the scatter for 
both families of E/SOs. 

In the mass-radius relation, blue-sequence E/SOs scat- 
ter toward larger radii than red-sequence E/SOs at 
fixed mass. Moreover, many red- and most blue- 
sequence E/SOs follow a shifted mass-radius relation be- 
low log M t /M ~ 9.7 (Fig. [9]), in the same regime where 
blue-sequence E/SOs become abundant. The trend to- 
ward larger radii for blue-sequence E/SOs is not dramatic 
in this regime — below M t , its K-S test significance is 
only ~2<j — but it is reproduced in both samples (Fig.[5]). 
To minimize the effect of stellar population differences on 
this result, we have adopted the if-band elliptical Kron 
radius from the 2MASS database as our size measure. 
This radius is typically ~2.5 x the elliptical half-light 
radius (and is in fact defined as 2.5 x the first moment 
of the light distribution), but in practice r^ ron displays 
less scatter than rff in the mass-radius relation. 

In the mass-cr relation, scatter increases dramati- 
cally below Mb (Fig. ITU]) , most notably for spirals and 
red-sequence E/SOs. Comparing residuals relative to 



the overall fit, K-S tests show that spirals and blue- 
sequence E/SOs are drawn from different parent popula- 
tions at >3cr confidence, whereas red- and blue-sequence 
E/SOs are consistent with the same parent population. 
Counting mid-sequence objects as blue- rather than red- 
sequence E/SOs slightly differentiates the two types of 
E/SO, yielding ~17% probability that both are drawn 
from the same population below Mb- To the extent that 
there are differences in this mass regime, blue-sequence 
E/SOs seem to have low stellar velocity dispersions more 
consistently, or equivalently, seem to lack scatter toward 
high dispersions. We suspect that the higher scatter in 
red-sequence a values points to greater dynamical diver- 
sity, especially since this scatter is not driven by errors: 
the four high-dispersion red-sequence E/SOs below Mf, all 
have very small errors on a. However, a larger sample 
would be needed to confirm the scatter difference statis- 
tically, and high-quality kinematic data sensitive to small 
a values are not available for the HyperLeda+ sample. 11 

The mass-cr relation of Fig.rTUlis not the Faber- Jackson 
relation between bulge mass and velocity dispersion, but 
rather the global relation between total stellar mass and 
central velocity dispersion. A galaxy's position in this 
global relation may depend on bulge-to-disk ratio and 
on how much the system is supported by rotation, in- 
cluding whether the bulge is dynamically cold, as ex- 
pected for a disky "pseudobulge " formed by disk gas in- 
flow and central star formation ([Kormendv fc Kennicuttl 
2004). For comparison, the right panel of Fig. ITOl shows 
the global relation between internal velocity and stellar 
mass, adopting the simplistic assumption that spirals are 
rotation-supported and E/SOs are dispersion-supported, 
and scaling the dispersions to equivalent circular veloc- 
ities by a factor of y/2 (e.g.. iBurstein et aT1 [l997). The 
four high-dispersion red-sequence E/SOs below Mb fol- 
low the relation defined by spiral galaxies and high-mass 
E/SOs, whereas the lower-dispersion E/SOs fall short in 
dynamical support, implying that they are probably sup- 
ported by_rotatioji(asjs commo this mass 
range: iDavies et al.lfl98l iBender et alill992ft . Unfortu- 
nately, in most cases the rotation velocities of the low- 
dispersion E/SOs are impossible to measure with our ex- 
isting data, due to abundant warps and asymmetries in 
their ionized-gas and stellar rotation curves as well as fre- 
quent radial truncation of the ionized gas emission lines 
(see also %WM- 

A tendency toward larger radii in blue-sequence E/SOs 
can be interpreted as consistent with either evolution 
from a more compact state via disk building or incom- 
plete evolution toward greater central concentration, as 
in a post-merger starburst phase. Lower cr's for blue- 
sequence E/SOs would be more naturally explained in the 
disk-building scenario, assuming that in the post-merger 
scenario the galaxy should have experienced dynamical 
heating (e.g.. iDasvra et al"1l2006f ). but a gas-rich merger 
could involve immediate inner disk regrowth via inflows. 

These results rely on stellar mass estimates that 
may vary systematically between subpopulations with 
different characteristic star formation histories, de- 

11 SDSS velocity dispersions are limited to logc > 1.8—1.9 due 
to instrumental resolution. Within slight differences in details 
of M, and a estima tion, our M*-cr relation agrees with that of 
Gallazzi ct al. (2006) above M b , but we find higher scatter in a at 
lower masses. 



Blue-Sequence E/SO Galaxies 



9 




log stellar mass (M SU(1 ) log stellar mass (M sun ) 

Fig. 5. — Subsamples of the NYU VAGC low-redshift sample with well-defined statistical properties. Contours in both panels are derived 
from the full sample (not shown), which spans redshifts 1000-15000 km s~ x and has an apparent magnitude limit of r = 18. (a) Apparent- 
magnitude limited VAGC subsample with r < 14.9. Light and dark gray points show red- and blue-sequence galaxies, respectively, (b) 
Volume-limited sample derived by including all VAGC galaxies down to M r = —17.15 with an upper redshift limit of 7500 km s _1 (light 
gray points). Dark gray points show fainter galaxies, demonstrating that this volume-limited sample is largely complete in stellar mass 
down to log M»/M ~ 9. 
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Fig. 6. — Comparison of stellar mass distributions for the two VAGC samples of Fig. [5] and the HyperLeda+ and NFGS samples. Red- 
and blue-sequence galaxies are shown in gray and black, respectively. For the two VAGC samples, open histograms show all selected 
galaxies, while filled/cross-hatched histograms show the subset having a 2MASS _fC-band flux with error <0.3 mag, matching the selection 
criteria used for the Hyper Leda+ and NFGS samples (§ I2.H . 



pending on the methods and assumptions used 
(jKannappan &: Gawiseii I2007D . Substituting stellar 
masses derived from the color- M^/L calibration of 
I Bell et all (|2003l ) yields greater scatter in the M*-a and 
M*-internal velocity relations and more strongly differ- 
entiates red- and blue-sequence E/SOs in these relations. 
The offset in radii between the two populations dimin- 
ishes under this substitution but does not disappear. 
Conversely, adoption of the stellar ma s ses de rived by 
Kannappan & Gawiser using iMarastonl (|2005h popula- 
tion synthesis models enhances the mass-radius offset 
while diminishing the mass-a offset. However, revised 
mass calculations using the more representative model 
grid of star formation histories described in § I2.3I yield 
very similar results for Maraston and Bruzual-Charlot 
models. 

4. CHARACTERISTICS OF BLUE-SEQUENCE 
E/SOS 

Here we offer an overview of blue-sequence E/S0 prop- 
erties, with special attention to properties relevant to 



evaluating the disk-building and fading-merger scenar- 
ios. 

4.L Morphologies 

Blue-sequence E/S0 morphologies depend strongly on 
mass (Fig. [3]). At high masses between Mb and M s 
(logM*/M Q = 10.5-11.2), many blue-sequence E/SOs 
resemble interacting galaxies or recent merger remnants. 
Of the 83 blue-sequence E/SOs in the HyperLeda+ sam- 
ple, we assigned seven the type "pec-e" to indicate 
spheroidal but otherwise unclassifiably disturbed mor- 
phology, and all seven have masses between Mf, and M s . 
The 16 highest-mass blue-sequence E/SOs include five 
with "pec-e" morphology and two more with prominent 
tidal features or dust lanes. The Mb~M s population is 
best probed by the HyperLeda+ sample because of its 
top-heavy luminosity distribution, but we find similar 
evidence for tidal features and interactions in two of the 
four highest-mass blue-sequence E/SOs in the NFGS sam- 
ple. At the same time, this mass range is also where the 
red and blue sequences converge, and some high-mass 
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Fig. 7. — Frequency of red- and blue-sequence E/SOs as a function of stellar mass in the HyperLeda-l- sample, (a) Raw numbers of 
objects, (b) Percentages of blue-sequence E/SOs among all E/SOs (solid line) and amo ng a ll galaxies (dotted line), using the NYU VAGC 
low-redshift catalog to completeness-correct the HyperLeda-l- sample as described in § 1 3 ■ 1 1 
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Fig. 8. — Distribution of HyperLeda+ and NFGS galaxies in 
normalized global environmental density vs. stellar mass. Black 
crosses mark blue-sequence E/SOs, gray dots mark red-sequence 
E/SOs, and small black dots mark late-type galaxies. Densities are 
expressed in units of the mean density of gal axies brighter tha n 
M B ~ -17 in the Updated Zwicky Catalog IIFalco et . al.|[l999h . 
smoothed on ~7 Mpc scales, using code adapted from N. Grogin 
(Grogin & Gcller 1998). In these units the densities of the Virgo 
and Coma clusters are ~4.9 and 7.4, respectively. To minimize 
edge effects, we plot only galaxies lying more than one smoothing 
length from the edge of the UZC and having redshifts cz > 1000 
and < 9500 km s with respect to the Local Group. 



blue-sequence E/SOs in the HyperLeda-l- sample appear 
quite normal; these may simply be red-sequence E/SOs 
affected by scatter in u — r. 

At lower masses where blue-sequence E/SOs are com- 
mon, strongly disturbed morphologies are rare. Most 
systems have fairly settled, disky morphologies. Faint 
tidal features and evidence of satellite interactions seem 
to be frequent but are hard to quantify (these features are 
not generally visible in the small, low-contrast cutouts 



in Fig. [3]) ; two quantitative but indirect measures of ex- 
ternal disturbance are discussed in § 14.31 In the NFGS 
sample, nearly all E/SOs with M* < M& have SO-SO/a 
types, regardless of sequence (Fig. [TJd), as most Es fall 
on the high-mass red sequence (c onsistent with the E 
and SO luminosity functions, e.g.. iBinggeli et aLl ll988'). 
In the HyperLeda-|- sample we find no unambiguous Es 
among 83 confirmed blue-sequence E/SOs (four galaxies 
are ambiguous between E and SO). Vetted types are not 
available for the red sequence for this sample. 

An analysis of the detailed morphologies of E/SOs is 
beyond the scope of this paper, but we note two possibly 
interesting facts that merit further investigation: (1) the 
ratio between SO and SO/a types is not obviously different 
between the red and blue sequences, based on the NFGS 
sample; and (2) bars, rings, and dust are evident in both 
red- and blue-sequence E/SOs (Fig.[3|). 

4.2. Colors, Gas, and Star Formation 

Intermediate to low mass blue-sequence E/SOs are 
bluer than red-sequence E/SOs in both their outer disks 
and their centers, with the strongest statistical difference 
being for outer-disk colors (measured between the 50% 
and 75% light radii). Fig. [TT] shows central and outer- 
disk colors for NFGS E/SOs, in the mass range below 
logM„/MQ = 10.7 where we have similar numbers of 
objects on both sequences. The two mid-sequence E/SOs 
are separated out. Comparing the central and outer- 
disk colors for this subsample, 7 of 14 blue-sequence 
E/SOs (50 ± 10%) have centers bluer than their outer 
disks, whereas only 2 of 13 red-sequence E/SOs are blue- 
centered (15~!l\ %). The disparity is greater if we limit 
the comparison to masses below M t , yielding percent- 
ages 861^% and 33+*4%, respectively. The two mid- 
sequence E/SOs are both blue-centered, with outer-disk 
colors intermediate between those of typical red- and 
blue-sequence E/SOs. Note that these blue central col- 
ors are not caused by AGN: the three known AGN in 
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Fig. 9. — X-band radius vs. stellar mass for (a) the HyperLeda+ sample and (b) the NFGS sample. The horizontal bar represents a 
typical 68% confidence interval in mass (i.e., from —la to +1<t). We use 2MASS Kron radii to emphasize underlying stellar structure; 
these radii lack catalogued uncertainties but do not seem to contribute much scatter. Lines in panel a represent least-squares bisector fits 
for E/SOs above and below Mt, illustrating the offset of blue-sequence E/SOs toward larger radii as well as the general shift in E/SO locus 
below the threshold mass. 
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Fig. 10. — Dynamical scaling relations, (a) Stellar mass vs. central velocity dispersion a for the NFGS sample. Velocity dispersions are 
measured within an r e /4 aperture. We exclude measurements with formal uncertainties >70 km s _1 or >cr/2 and add 5% in quadrature 
(typical repeatability) to the formal uncertainties to compute the error bars shown. The line is a least-squares bisector fit to all the data. 
(b) Stellar mass vs. internal velocity for the NFGS sample. E/SOs are plotted with velocity dispersions as in panel a, here multiplied by 
\/2 to scale with rotation velocities (e.g.jBurstcin ct al. 1997). Late-type galaxies are plotted with rotation velocities, equal to half the 
inclination-corrected Wv pmm measure of !Kaimappaii^t*EdTl2002l) . Strongly asymmetric or truncated rotation curves are excluded by the 
criteria of Kannappan & Barton (2004). Lines mark key mass and velocity scales, as indicated. 



our sample are all in red-centered systems (one Seyfert 
2 and tw o LINERs, identified with NFGS nuclear spec- 
troscopy: Uansen et al.1 l2000a|) . The high frequency of 
blue-centered profiles is intriguing as possible evidence 
for interaction-triggered central star formation (§ 14.3ft . 
Nonetheless, all of these systems have blue outer disks 
as well, as also seen in blue-sequence E/SOs that are not 
blue-centered. 

Fig. [12] displays gas detections and gas-to-stellar mass 
ratios for E/SOs in the NFGS sample. Stars mark five 
systems with counterrotating or polar-ring gas, to be an- 
alyzed in detail in § [SJ Although the archival HI data are 
incomplete and the ionized gas detections are somewhat 



biased against red-sequence E/SOs, 12 one can still see a 
pattern of more frequent gas detections below M t on the 
red sequence and below Mb on the blue sequence, re- 
flecting an overall i ncrease in gas fractions for the entire 
galax y population ([Kann appan 2004; Kannappan fc Weil 
l2008h . 

At a given stellar mass, blue-sequence E/SOs have more 
gas than red-sequence E/SOs, with atomic-gas-to-stellar 
mass ratios ranging from ~0.1 to >1. However, the dis- 
tributions overlap a bit and do not obviously preclude 
that some red-sequence and mid-sequence E/SOs may 

12 Spectra for these systems were obtained in a wavelength range 
around H/3 rather than the stronger Ha line. 
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evolve toward the blue sequence, either gradually or in 
intermittent bursts, especially if their environments are 
favorable to fresh gas accretion. 

To assess the potential for morphological transforma- 
tion in blue-sequence E/SOs, via either disk building or 
central mass growth, we consider several measures of star 
formation summarized in Fig. 1131 We compute the spe- 
cific star formation rate (SSFR), expressed in units of 
percentage growth per Gyr at the current star forma- 
tion rate, as well as two star formation timescales: the 
stellar mass doubling time, equal to current stellar mass 
divided by current SFR, and the atomic gas consump- 
tion time, equal to atomic gas mass divided by current 
SFR. These are approximate, instantaneous measures, 
uncorrected for any future decline in SFR, recycling, or 
infall. Star formation rates are derived from integrated 
Ha spectral line data from the NFGS database, origi- 
nally obtained by s canning the slit across each galaxy 
(|Jansen et al.ll2000af) . We correct for extinction and con- 
vert to t he common scale cal ibrated against IRAS-based 
SFRs bv lKewlev et all ((2002). Where available, we adopt 
the catalogued "Ha c<lrr " SFRs provided by Kewley et al., 
and we compute SFRs for additional galaxies in our sam- 
ple using integrated Ha fluxes provided by R. A. Jansen 
(private communication, 2001). All SFRs are converted 
to our H a and IMF. 

From these data, blue-sequence E/SOs have SSFRs 
comparable to those of spiral galaxies (Fig. [T3k). Rates 
of transformation range from ~2-20% per Gyr. While 
SSFRs generally increase toward lower mass for both 
spirals and E/SOs, the two most rapidly transforming 
E/SOs have masses between M t and Mj,. The highest- 
mass systems appear quiescent, as do most red-sequence 
E/SOs, but one mid-sequence system is evolving at a rate 
of nearly 15% per Gyr. 

Most blue-sequence E/SOs, as well as two red-/mid- 
sequence E/SOs, have enough gas to form stars for several 



Gyr. Fig. 113b shows doubling time vs. gas consumption 
time. To the left of the line of equality, galaxies can- 
not double their stellar masses without fresh gas accre- 
tion (unless the unmeasured molecular gas mass is large) . 
Four blue-sequence E/SOs and one mid-sequence E/SO 
combine doubling times less than a Hubble time with 
gas consumption times that are either longer or <2-3x 
shorter, implying significant evolution even in an unreal- 
istic closed-box scenario. Notably, the two-blue sequence 
E/SOs whose gas consumption times match or exceed 
their stellar mass doubling times are both engaged in 
strong interactions with smaller, but substantial, com- 
panions. 

While the above analysis suffers from incomplete data 
(in particular, one of the two high SSFR systems in 
Fig. [T3k lacks HI data in in Fig. H3b). our results con- 
firm the plausibility of the basic premise of morpholog- 
ical evolution. Further analysis of potential structural 
changes will require resolved gas and star formation data, 
which will soon be available from Spitzer, GALEX, GBT, 
CARMA, and VLA programs underway (Wei et al., in 
prep., and Moffett et al., in prep.). 

4.3. Interaction Status 

We have already established the importance of strong 
interactions for high-mass blue-sequence E/SOs, based on 
their morphologies (§ I4.1[) . However, weaker interactions 
such as satellite accretion events are not easy to identify 
cleanly based on morphology. Here we offer two other 
lines of argument to support our subjective impression 
that minor interactions are important for blue-sequence 
E/SOs below M s . 

First, the close correspondence between blue-sequence 
E/SOs and blue-centered color gradients found in § 14.21 
may reflect star formation triggered by interactions. 
In an earlier study o f NFGS galaxies of all types, 
iKannappan et al. (2004) showed that blue-centered color 
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gradients correlate with morphological peculiarities sug- 
gestive of minor mergers and argued that tidal forces in 
these accretion events can trigger disk-gas inflows and 
central star formation, p ossibly leading to the formatio n 
of disky "pseudobulges" (jKormendv fc Kennicuttl feutM). 
Most of our blue-centered E/SOs were not flagged in this 
prior study as morphologically peculiar, perhaps because 
by definition E/SOs are fairly regular. Therefore if blue- 
centered color gradients in E/SOs reflect recent interac- 
tions, the disturbances are likely mild. 
As evidence that such mild disturbances may in- 



deed be at work, we find that below Mft, 6 of 11 
(55 ± 11%) blue-sequence E/SOs with ionized gas ro- 
tation curves in the NFGS database show strong rota- 
tion curve asymmetries (>10% in the asymmetry met- 
ric of iKannappan fc Barton! 2004). Five of these six are 
also blue-centered, although that correspondence is not 
statistically significant given the small numbers. More 
significantly, the rate of strong asymmetries for blue- 
sequence E/SOs is >3 times that seen for the general 
galaxy population in this mass range (16 ± 2%). Unfor- 
tunately, too few red-sequence and mid-sequence E/SOs 
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in our sample have the extended emission lines necessary 
to measure rotation curve asymmetries: only 5 have gas 
rotation curves, of which 3 are truncated at <0.9r e . 

Of course, companion statistics would be valuable 
in assessing the role of external disturbances. While 
small satellite data are not yet available, we have iden- 
tified large companions within 300 kpc and 300 km s _1 
of all NFGS galaxies that lie inside the boundaries o f 
the Updated Zwicky Catalog fUZC; llMcoetal] [l999h . 
These data reveal similar large-companion rates 13 among 
red-sequence E/SOs, blue-sequence E/SOs, and late-type 
galaxies, subject to data incompleteness. If anything, the 
rate of large companions may be marginally lower among 
blue-sequence E/SOs. Clearly this question would benefit 
from further investigation. 

5. COUNTERROTATING AND POLAR DISKS: 
SECONDARY STELLAR DISK GROWTH 

An especially distinctive feature of blue-sequence 
E/SOs is their association with counterrotating gas and 
polar rings. As shown in Fig. [TJ all five NFGS E/SOs 
with decoupled gas lie on or near the blue sequence. 
Whether and how these galaxies may form stars in their 
decoupled gas is highly relevant to understanding the 
evolutionary destiny of blue-sequence E/SOs. Here we ex- 
amine these case studies in some detail, deferring a more 
general discussion of formation scenarios to Section [5J 

5.1. The Polar Ring: Triggered Stellar Disk Growth 

The polar ring galaxy UGC 9562 is the bluest and 
lowest-mass of the five decoupled gas systems and con- 
tains an obvious secondary stellar disk, i.e., the ring. 
This system is shown in Fig. [3] with (logM*, color) near 
(8.9, 1.3). UGC 9562 has rec ently interacted with a com- 
panion, but lCoxeT al. (2001) argue based on stellar pop- 
ulation analysis and HI gas morphology that the interac- 
tion has simply triggered star formation in a preexisting 
gas ring. Fig. [14] shows new evidence for t he preexist- 
ing ga s scenario, based on data reported in IGuie et aLI 
(2005). We plot the ionized gas and stellar kinemat- 
ics for UGC 9562 and its companion overlaid on the HI 
position- velocity cut along the ma jor axis of the g as be- 
tween the galaxies (see Fig. 3 of tCox et al.l l200i[ ). As 
observed by Cox et al., the HI connecting the galax- 
ies appears to form a distinct structure from the ring. 
Moreover, our new data show that near the center of 
UGC 9562, its ionized gas has higher line-of-sight veloc- 
ity than its stars, while the companion's ionized gas and 
stars both have somewhat lower line-of-sight velocities 
than UGC 9562 's stars. These observations are not con- 
clusive but most naturally suggest a recent flyby event in 
which previously accreted gas was tugged and disturbed, 
rather than an event in which the lower-velocity system 
directly donated the high-velocity gas. 

These results support the preexisting gas scenario, 
which implies that (a) a new stellar disk is forming 
in situ in gas that has accreted onto a previously 
formed E/S0 and that currently dominates the total 
mass (Mh, i+Uc/M* ~ 3), and (b) the mechanisms and 

13 Here "large" signifies "as much as ~1 mag fainter," because 
the UZC was derived from the CfA 2 redshift survey catalog, whose 
limiting magnitude is 1 mag fainter than the CfA 1 survey from 
which the NFGS was selected. 



timescales of gas and stellar disk formation may be dis- 
tinct, with stellar disk growth being potentially bursty 
and possibly requiring a trigger. Cosmological simula- 
tions suggest a picture in whic h massive polar ring s may 
form via cold- mode accretion (jMaccio et al.l[20 0lf). con- 
sistent with observations that favor a slow gas accretion 
scena rio for another polar ring, NGC 4650A ( Iodice et al.l 
[20061) . 

5.2. The Counterrotators: Past and Present Disk 
Growth 

The other four E/SOs marked with stars in Fig. [T2l were 
identified bvlKannappa n fc Fabri cant (20Ql|) in a system- 
atic search for NFGS galaxies with counterrotating gas 
and stars. Fig. [T5l shows new, high S/N stellar absorption 
line cross-correlation profiles for these galaxies, revealing 
their stellar velocity substructure as a function of slit po- 
sition. These data were obtained in May 2001 with the 
Blue Channel spectrograph on the MMT, using a config- 
uration with 1'.'2 binning and a ~ 40 kms" 1 resolution 
centered on the Mg I triplet at 5175A. The spectra are 
identical to those in the standard NFGS database, ex- 
cept that deep ~1 hr exposures enable us to search for 
secondary stellar velocity components. We measure the 
cross-correlation profile as a function of velocity at each 
slit position using xcsao in the rvsao package of IRAF 
([Kurtz et al.l fl9 9Z) . correlating against non-rotating G 
and K giant stars observed in the same instrumental con- 
figuration. Here we consider what may be learned from 
the raw cross-correlation profiles, deferring full modeling 
of the line-of-sight velocity distributions. 

Secondary stellar components that rotate with the gas 
are clear in two systems: (1) the very disky elliptical 
NGC 7360, which shows resolved X-structure in its cross- 
correlation profiles, indicating oppositely rotating stellar 
components; and (2) the low-mass SO NGC 3011, which 
shows stars and gas that rotate together in the central 
region, but in opposite senses at large radii, with profile 
asymmetries suggesting that the stellar component that 
rotates with the gas may be present but subdominant 
and unresolved at larger radii. Profile asymmetries in a 
third system, NGC 5173, likewise suggest the possibility 
of an extended, unresolved secondary stellar component 
rotating with the gas. While the asymmetries are not 
conclusive by themselves, a young stellar disk is detected 
morphologically in this galaxy, whose primary round E 
component hosts an inclined large- scale spiral disk with 
blue knots ([Vader fc Vigrouxlll99lH . 

In contrast with these three systems, the fourth coun- 
terrotator, UGC 6570, shows no hint of a secondary stel- 
lar component, although its outer velocities seem to re- 
flect a strong warp or distortion, perhaps from a re- 
cent interaction. This galaxy may be the exception that 
proves the rule regarding secondary stellar disk forma- 
tion. We suspect that its gas arrived too recently to have 
had time to form stars, because spectra at different posi- 
tion angles indicate tha t the gas rotates in a tilted plan e 
relative to the stars (Kannap pan fc Fabrica nt 2001), 
which is an unstable and short-lived configuration. If 
such an accretion event only recently pushed UGC 6570 
off the red sequence, then its status as a mid-sequence 
E/S0 between sequences may change as star forma- 
tion proceeds, so that it may soon join the other po- 
lar/counterrotating gas systems on the blue sequence. 



Blue-Sequence E/SO Galaxies 



15 



1350 - 



1300 



I 1250 - 



1200 



1150 



222.70 



1350 FT^ 



UGC 9562 



UGC 9560 



1240 



1200 



1180 . gas 
□ stars 



1250 - 



1200 - 



nil 



' [] 



ring/main gas • » 
ring/main stars □ A 



-20 -10 10 
arcsec along slit (NOT TO SCALE) 



-20 -10 10 20 30 
arcsec along slit (NOT TO SCALE) 



222.72 222.74 222.76 222.78 

RA (degrees, 2000) 



222.80 



222.82 



222.84 



Fig. 14. — Schemat ic co mparison of velocities in different components of the companion galaxies UGC 95 60 and UG C 9562 (the polar-ring 
galaxy discussed in § 15.11 1. Grayscale approximates the lowest surface brightness HI contour detected bv lCox et al.l (2001). Open triangles 
and squares mark stellar absorption-line velocities, while filled triangles and circles mark ionized gas emission-line velocities from H/3 and 
the nearby [OIII] lines. Inset panels are positioned correctly in position and velocity, but the spatial scales are expanded to show detail. 



It is empirically clear that three of our four counterro- 
tators have managed to form significant secondary stel- 
lar components, at least one of which seems to be gain- 
ing dominance from the inside out (NGC3011), and an- 
other one of which has unmistakably substantial mass 
(NGC7360). Yet these systems have only ~10-20% gas 
fractions and among the lowest SSFRs in Fig. [131 Bulk 
accretion of both gas and stars from an outside galaxy 
into a disk-like secondary component seems unlikely to 
explain these galaxies, considering that the mid-sequence 
system UGC 6570 appears to contain externally accreted 
gas without accompanying stars (as we also infer for the 
polar ring system in § 15. 1[) . Also, in most of these sys- 
tems the counterrotating gas is extended, sometimes well 
beyond the stars (NGC 7360 shows HI and Ha emission 
out to >3x its optical radius; Kannappan, Matthews, 
Christlein, et al., in prep). We suspect that the dis- 
parity between past and current star formation reflects 
a history of intermittent bursts of enhanced stellar disk 
growth, occurring independent of gas disk growth. Also, 
the presence of coextensive counterrotating stellar popu- 
lations may exacerbate such intermittency, if the galaxy 
is stabilized against forming spiral arms or other struc- 
tures, so that only strong interactions can revive star 
formation (Kannappan, Matthews, Christlein, et al., in 
prep). 

Our detection rate of 2 definite + 1 probable counter- 
rotating stellar disks in a sample of 4 gas-stellar coun- 
terrotators m a y see m to contradict the null result of 
iKuiiken et al.l (|1996f) . who surveyed 28 SO galaxies for 
counterrotating stellar disks and found none. However, 
we note that they analyzed only 4 systems with gas- 



stellar counterrotation, so the samples that are relevant 
to compare are tiny and may differ in crucial param- 
eters such as mass and color. The fact that Kuijken 
et al. could not find any counterrotating stellar disks in 
SOs without gas might also have another explanation: if 
counterrotating stars reflect disk regrowth, then SOs with 
counterrotating stars may be on their way to becoming 
spiral galaxies. If this conversion is efficient, the condi- 
tion of still having SO morphology but having built up a 
detectable counterrotating stellar disk may be an unusual 
state, becoming long-lived only when the disk regrowth 
process is quenched prematurely. 

6. DISCUSSION 

While individual blue-sequence E/SOs have been stud- 
ied for years without being identified as a class, recogniz- 
ing these galaxies as a coherent population opens the way 
to understanding their role in the story of galaxy evolu- 
tion. Below, we break down blue-sequence E/S0 forma- 
tion scenarios as a function of mass, then take a closer 
look at the disk regrowth scenario in light of other stud- 
ies. We close by considering how these results may affect 
our understanding of galaxy evolution from z = 1 — 0. 

6.1. Evolution in Three Mass Ranges 

Taken together, our results suggest that blue-sequence 
E/SOs are a composite population, reflecting different 
physical processes in different mass regimes. We distin- 
guish three regimes: above the shutdown mass M s , below 
the threshold mass M t , and in the intermediate range 
centered on M&. The intermediate regime is the primary 
mass regime of large spiral galaxies and as such defines 
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a coherent population; however, blue-sequence E/SOs in 
the upper and lower halves of this regime show distinct 
evolutionary trajectories. 

Above M s , blue-sequence E/SOs are virtually nonex- 
istent, as the blue sequence itself largely disappears. 
From M b up to M s (log M«/M - 10.5-11.2), there is 
a clear decline in the frequency of blue-sequence E/SOs 
(Fig. [7|), and their morphologies often show signs of 
strong disturbance. Their <1% abundance relative to 
other massive galaxies is roughly comparable to the fre- 
quency of close pairs, suggesting an approximate corre- 
spondence to the merger rate, depending on relative time 
spent in the close-pair and fading-starburst p hases (e.g., 
iDe Propris et alj|2007t rPomingue et al.ll2009h . These re- 
sults point to major and substantial minor mergers as the 
probable source of most high-mass blue-sequence E/SOs. 
Because cold gas is generally scarce in E/SOs with masses 
above Mb — both within the galaxies and in their typi- 
cally dense environments (Figs. [8] and [T2|) — high-mass 
blue-sequence E/SOs are likely to exhaust all available gas 
in their current post-merger phase and migrate perma- 
nently to the red sequence ("quenching mergers"). How- 
ever, limited disk building is also occurring in this mass 



regime (e.g., NGC7360, § [5j) . 

At the other extreme, below the threshold mass 
log Mt/M ~ 9.7, blue-sequence E/SOs show little re- 
semblance to canonical major-merger remnants. Here 
blue-sequence E/SOs make up 20-30% of all E/SOs and 
~5% of all galaxies. The coincidence of their emergence 
with a gene ral rise in gas richness for the entire galaxy 
population (|Kannappanll2004l iKannappan fc Weiil200ct ) 
strongly suggests that the availability of cold gas is fun- 
damental to their existence. 

Intermediate-mass E/SOs from down to M t are of- 
ten gas rich as well, at least on the blue sequence, and 
our two most rapidly growing E/SOs lie in this mass range 
(Fig. [T3|) . Available environmental data reveal a strong 
shift toward lower density environments for E/SOs be- 
low Mb, down to at least 10 9 Mq, so fresh gas accretion 
seems plausible (Fig. [8] recall that cluster dEs become 
abundant at lower masses). The shift in environments 
applies to both red- and blue-sequence E/SOs, and the 
two populations roughly track each other in numbers be- 
low Mb (Fig. [7] and § 13. 1[) , but red-sequence E/SOs from 
Mb down to M t can be gas poor nonetheless. In this mass 
range, red-sequence E/SOs seem to show greater diversity 
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in V/a than blue-sequence E/SOs (within small number 
statistics; see Figs. [10]). A shifted M* -radius relation be- 
comes dominant below M t for E/SOs on both sequences, 
but some red-sequence E/SOs fall on an extension of the 
higher-mass M*-radius relation (Figs. [9]). Blue-sequence 
E/SOs make up 5-7% of all E/SOs and nearly 2% of all 
galaxies between Mb and M t , potentially reflecting evo- 
lution in both directions between the sequences rather 
than uni-directional quenching. 

6.2. Forming and Re- Forming Disks 

Our data suggest that blue-sequence E/SOs in mass 
and environment regimes with abundant cold gas are 
most likely engaged in building disks. Blue outer disks 
are the strongest statistical signature of blue-sequence 
E/SOs (§ 14. 21) . and a number of them show evidence for 
secondary stellar disk growth in counterrotating or polar 
gas (§©. At the same time, the fact that blue-sequence 
E/SOs are more compact and dynamically hotter than 
spirals, and more similar to red-sequence E/SOs than 
spirals in scaling relations, suggests that most of these 
E/SOs have passed through a violent merger phase at 
some point in their formation history (§ 13. 2|) . Other 
processes that could yield compact structure, e.g., ha- 
rassment, are less prevalent in the modest-density en- 
vironments seen for the sub-Mf, E/SOs in our samples. 
Moreover, these processes would not be expected to pro- 
duce blue outer disks or secondary disk growth. Thus the 
extended blue disks seen in blue-sequence E/SOs plausi- 
bly represent true disk regrowth over merger-formed early 
type galaxies, as envisioned in hierarchical models (e.g., 
ISteinmetz fc Navarrdl2002l; iGovernato et al.ll2007t) . 

A variety of hints in the literature point to the 
likelihood of disk regrowth around E/SOs, especially 
SO galaxies. Many SOs contain gas that appears 
to have an exte rnal origin: HI rings ti l ted o ut of 
the stella r disk (Ivan Driel fe van Woerdenl I1991I ). po- 
lar rings dBettoni et alJ l2001h . or counterrota ti ng ga s 
(jBertola et all Il992f) : see also ISage fc Welchl (|2006h . 
These extended HI disks and rings are often associ- 
ated wit h star formation and/or faint extended stel- 
lar disks (lHawarden et alJll98lHvan Gorkom et alJll987t 
iNoordermeer et al.l 120051: 1 Jeong et al.ll2007D . "Antitrun- 
cated" disks with a low surface brightness component 
that dominates at large radii are also reported in SOs 
(|Erwin et alJ [2005). Recent star formation in E/SOs 
may be substantial, involving >10 % of the stellar mass , 
especially in field SO galaxies (lAnnibali et al.l 120071: 
Nolan et aT]|2007tlKavirai et alj|2007bt ISchawinski et all 
2007T see also lGallazzi et al. I 20061: iDonas et al.ll2007l )7Tn 
polar ring galaxies and other systems where secondary 
disks can be isolated, these co mponents can be very mas- 
sive (§ EJ IHtonLitSllMJ). 

We speculate that disk-building E/SOs may (re)join 
the spiral sequence when secondary disks become domi- 
nant, perhaps helping to explain the surprising scarcity 
of SOs with two c ounterrotating stellar disks (§ 15.21 
iKuijken et al.l Il996h and the existence of peculiar ob- 
j ects like the prolate-bulg e Sbc galaxy UGC 10043 
([Matthews fc de Griisl [2004). which could be a former 
polar ring galaxy whose ring has become a large disk 
around the central E/S0. To grow into spiral galaxies, 
blue-sequence E/SOs must grow on average by a factor 
of —2 in radius (Fig. [9]), representing a huge increase in 



mass. Thus (re)forming a late-type system on a cosmo- 
logically relevant timescale will require two things: abun- 
dant gas and efficient star formation. 

Efficient star formation does n ot seem to characteriz e 
the most massive gas- rich SOs (jOosterloo et al.l l2007t ). 
which include galaxies suc h as Malin 1 (recently ex- 
posed as an SO at heart bv iBarthl [200l and NGC3108 
(|Hau et aT1l2008f ). Inflow mechanisms appear to be in- 
effective for E/SO s with dispersions a ~ 230 km s _1 
(|Serra et al.ll2008l ). equivalent to M* > M s (Fig. ITU)) , 
Nonetheless, minor mergers may help to drive disk gas 
inward in i ntermediate-mass E/SOs, fueling residual sta r 
formation (Kauffm ann et aill2007l : lKaviraj et al.l l2007a). 

The real action, however, starts below Mb and in- 
creases below M t , down to <10 9 Mq. Here blue-sequence 
E/SOs are likely to have formed in mergers of very gas- 
rich pro genitors (i.e., galaxies typical of this mass r egime: 
Fig.lT2l lKannappadl2004HKannappan fe Weill2008h . The 
remnants would have inh erited substantial gas through 
the fallback of tidal debris ([Barne s 2002; Robe rtson et alJ 
120061: IStewart et all 120091) . Moreover, the sub-M t mass 
regime may permit continued growth by e f ficient cold- 
mode gas accretion dBirnboim fc Dckcl 2003; Keres et all 
120051: iDekel fc Birnboimll2006f ). Subsequent interactions 
can also supply fuel for star formation, via direct accre- 
tion of gas-rich satellites, gas transfers, or compression 
of ambient gas. 

Star formation becomes remarkably efficient for E/SOs 
below Mb, and the two most strongly evolving blue- 
sequence E/SOs in the NFGS, whose current SFRs cor- 
respond to —20% mass growth per Gyr, have Af* ~ 
1O 1O M0 (Fig.[T3|. This result may seem counterintuitive, 
since below M t , most disk galaxies are late- type dwarfs 
(Sd, Sm, Im), which have notoriously inefficient star for- 
mation compared to classical large spirals. However, re- 
cent work suggests that this inefficiency may be linked 
to low surface mass dens i ties ra ther than to low masses 
per se. iKauffmann et al . ( 2006) report that evidence of 
higher-intensity, shorter starbursts increases with stel- 
lar surface mass density /i* for disk-dominated galaxies, 
reaching a peak at a characteristic /i* — 3 x 10 s Mq 
kpc -2 before falling off for higher-mass/more spheroid- 
dominated systems as the star formation rate per unit 
stellar mass <SFR>/M* decreases. The characteristic 
/i* corresponds to a characteristic concentration index 
C r ~ 2.5-2.6 (where C r is defined by the ratio of the 
radii containing 90% and 50% of the light in the SDSS 
r band). Galaxies near the characteristic /i* and C' r are 
common from M t up to just below M s (see Kauffmann 
et al.'s Figures 1 and 2). 

Intriguingly, blue-sequence E/SOs have nearly opti- 
mal C r : based on SDSS data for 73(50) NFGS galax- 
ies below Mb(M t ), blue-sequence E/SOs have mean C r — 
2.74(2.54), with standard deviation —0.45, whereas late- 
type galaxies have C r — 2.31(2.28) with standard devia- 
tion -0.40, and red-sequence E/SOs have C r = 3.01(2.86) 
with standard deviation —0.30. Kauffmann et al. inter- 
pret the peak in burst intensities near the characteristic 
/!* in terms of short gas consumption times for accreted 
clouds (equivalent to high star formation efficiencies), 
because they see fairly constant <SFR>/M* over the 
range of /x* below the characteristic value. However, con- 
stant <SFR>/M* implies relatively lower <SFR>/M tot 
for gas-rich late- type galaxies, a non-trivial distinc- 
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tion in the sub-M t re gime (Fig. [H iKannappanl [2004; 
iKannappan fc W ei 2008). Therefore, the Kauffmann et 
al. results are consistent with both higher <SFR>/Mt ot 
and higher-intensity bursts for blue-sequence E/SOs com- 
pared to spiral or red-sequence E/SO galaxies, apparently 
as a direct corollary of having surface mass densities near 
the characteristic value. 14 

It would seem that most blue-sequence E/SOs (and 
likely some red-sequence E/SOs and non-E/SO blue com- 
pact galaxies) lie in a sweet spot for disk building, with 
surface mass densities ideal for efficient star formation 
and masses low enough for abundant gas accretion. We 
conjecture that passing through a blue-sequence E/SO or 
similar phase in this sweet spot may enable a galaxy to 
escape the low surface density dwarf regime and form a 
large spiral disk. This conjecture is consistent with obser- 
vations that massive disk galaxies, no matter how thin or 
low surf ace brightness, always seem to show a bulge com - 
ponent (|Dalcanton et all [20041 : iSpravberrv et al.lll995fl . 
The blue-sequence E/SOs we have analyzed are a natu- 
ral progenitor population for large spiral galaxies with 
masses up to ~-M&, and this evolutionary link is likely 
to extend to higher masses at higher z (see § 16.3(1 . How- 
ever, the physics of this transformation, and particularly 
the growth of the disk beyond the radius of the original 
E/SO, is not yet clear. 

As the bulges of many late-type galaxies are disky, 
any late-type regeneration process should be able to 
form disky bulges. We may see this happening in 
blue-sequence E/SOs, whose frequent blue-centered color 
gradients and kinematic disturbances plausibly reflect 
episodes of pseudobulg e growth triggered by we ak inter- 
actions (§ 14.31 see also IKannappan et al.ll2004T ). Exter- 
nally or internally driven gas inflows can gradually build 
a larger pseudobulge over a classical bulge, as is seen in 
the d etailed HST decomposition of two SOs (|Erwin et alj 
|2003[ ) . In fact, pseudobulges are quite common in SOs 
(e.g jLaurikainen et ai1l2006f) . and iBarwav et alj (|2007f ) 
find evidence for a sharp transition in the scaling between 
SO bulge and disk radii below Mr- ~ —24.5 (in our sys- 
tem, equivalent to logM^/M© ~ 10.8-10.9), which they 
interpret as evidence of pseudobulge growth. We sus- 
pect that nearly all blue-sequence SOs have pseudobulges, 
given their modest velocity dispersions, whereas the large 
scatter in red-sequence E/SO velocity dispersions proba- 
bly indicates a variety of bulge types. Only systems in 
which the central spheroid is protected from gas inflows, 
as in a polar ring configuration, might preserve the orig- 
inal spheroid while building just the outer disk (as for 
the aforementioned prolate-bulge galaxy UGC 10043). 

The timescale for disk regrowth is probably longer 
than a typical merger timescale, given that most blue- 
sequence E/SOs show relatively regular morphology, 
and that their frequency in the galaxy population ex- 
ceeds that of close pairs (except near M s ; Fig. [7] 
and § 16. ip . However, small bursts of growth must 



14 Exactly what sets the optimal surface density is unclear in the 
absence of information on gas content, fi ga s, or the timing of gas 
accretion relative to starbursts. The population-averaged trends 
seen by Kauffmann et al. must include a variety of physics, as em- 
phasized by the falloff in burst intensities above the characteristic 
fi„. Both local disk dynamics and global effects of gas-rich minor 
or major mergers may link higher star formation efficiencies with 
higher fi t , while quenching mergers may work the other way. 



be common, as >50% of blue-sequence E/SOs show 
blue-centered color gradients and/or kinematic asym- 
metries (§ 14.31 14.21 and [5]). Bursty disk growth ap- 
pears to be the norm for large disk galaxies, if the 
star formation and accretion histo ries of the Milky 
Way and M31 are representative (Rocha-Pint o et al.l 
20001 Ide la Fuente Marcos fc de la Fuente Marcos! 12004} 



Helmi et al.l 12006 1 : llbata et all 120051: lElmegreen fc Scald 
20061). In fact. lHammer et all (j2007l) argue that the 
Milky Way's history may b e considered unusually quiet 
(but see Uames et al.l |2008|) . A variety of evidence sug- 
gests that disk galaxies receive regular fresh gas i nfall 
sufficient to fuel episodic bursts (|Sancisi et al.ll2008h . 

The likelihood of intermittent growth and/or dust- 
enshrouded star formation implies that some former or 
future blue-sequence E/SOs should fall on the red se- 
quence. Consistent with this picture, we find regimes 
of overlap in red-and blue-sequence E/SO properties (gas 
fractions, radii, velocity dispersions), with substantial di- 
versity among red-sequence E/SOs. With existing data, 
we cannot measure rotation velocities for the dynami- 
cally cold red- and blue-sequence E/SOs in our sample 
(§ 13. 2p . but many have stellar masses sufficient to join 
the spiral-galaxy Tully-Fisher relation above M t (i.e., 
above V ro t ~ 120 km s^ 1 ), provided they grow large 
rotating disks (Fig. [TOJ) . One such blue-sequence SO, 
UGC 12265N, has enough gas to double its stellar mass 
and is engaged in a strong interaction with a smaller 
but substantial companion, driving stellar mass growth 
at a rate of ^20% per Gyr (§ I4.2p . Because the even- 
tual merger should b e too gas rich to destroy the disk 
(Hopk ins et al.ll2009f ). this galaxy will almost certainly 
develop spiral structure as the disk forms stars and re- 
sponds to tidal forces in the interaction. Such an event 
may mark the birth of a large spiral galaxy. 

6.3. Implications for Galaxy Evolution from z = 1 — 

Most of the above discussion refers to a snapshot at 
z = 0. Given the downsizing of the bimodality mass, 
e.g., as seen in the downsizing of the crossover mass for 
E/SO and late-type ma ss functions from c lose to M s at 
z = 1 to today (B undv et al.1 120051 ) . we suspect 

that the gas-richness threshold mass M t has also evolved 
from z = 1 — 0. We therefore expect z ~ 1 analogues 
of blue-sequence E / SOs to show strong disk-building ac- 
tivity up to the equivalent threshold mass at z ~ 1, 
which might naively be expected to fall <1 dex below 
the z ~ 1 crossover mass, i.e. near a few x 10 10 M© 
(similar to today's Mb). Conversely, a noticeable frac- 
tion of z = galaxies with masses up to M s may repre- 
sent the end states of high-redshift blue-sequence E/SOs 
(or dusty red-sequence analogues: lHammer et al1l2005f) 
that have rebu i lt disk s. Consistent with this expectation, 
ISargent et ail (|2007l ) report from the COSMOS survey 
that the number density of large disks with bulges has 
doubled since z = 1, whereas they see no increase, and 
possibly a small decrease, in the number density of bulge- 
less disks at z = 0. These results point to the importance 
of forming bulges before building large disks, i.e., merg- 
ing gas-rich bulgeless late types to form E/SOs may be 
an enabling step in large disk galaxy formation (as also 
suggested by the star formation efficiency considerations 
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with masses up to a higher threshold mass perhaps near a 
few x 10 10 Mq, must certainly exist in high-z samples of 
blue spheroids, blue-co red spheroids, and luminous blue 
compact galaxies (e.g. Jim et alj|2001t IMenanteau et al 
20011: iTreu et all 120051 : lElmegreen et all 120051 : iLee et al. 
2006t ). Type-dependent luminosity function analyses 
reveal abundant blue E/SOs, whose luminosity distri- 
bution parallels that of late- type Irr /Pec systems, the 
analogues of today's d war f galaxies ilCross et alj l2004t 
IMenanteau et all |2006() . iKavirai et all (|2008l ) report 
a split in the color distribution of z ~ 0.5-1 E/SOs, 
with a "minor but significant peak" on the blue se- 
quence. Several authors have suggested links between 
blue co mpact systems and the bulges of today's large 
spira l s jKobulnicky fc Zaritskvl 119991 : iBarton fc van Zeel 
120011: iPuech et alj I2OO60 . However, both blue/blue- 
cored spheroids and blue compact galaxies are compos- 
ite populations, containing galaxies representing a wide 
r ange of masses , environments, and formation scenario s 
(|Lee et al.ll2006l; iNoeske et afll2006t IBarton et afll200ll . 
Even at z = 0, we have seen that evolutionary scenarios 
for blue-sequence E/SOs are a sensitive function of mass. 
Also, the bulk of the systems we have described reside 
in field environments, whereas blue or blue-cored E/SOs 
in clusters are more likely to reflect transient states as a 
prelude t o quenching by harassment, stripp i ng, or stran- 
gulat ion (|Koo et all 119971 : iRose et alJfeOOli iLisker et all 
l2006t) . Moreover, the abundance of AGNs among blu e 
spheroids is disputed (|Lee et al.ll200llTreu et al.ll2005f K 
Quantitative assessment of the disk regrowth picture will 
therefore require careful decomposition of high-z samples 
to isolate analogues of the disk-building blue-sequence 
E/SOs analyzed here. 

7. CONCLUSION 

We have examined the existence, properties, and evo- 
lutionary trajectory of a substantial population of mor- 
phologically defined E/SOs that fall on the blue sequence 
in color-stellar mass space, i.e., the usual locus of late- 
type disk galaxies. Our analysis of blue-sequence E/SOs 
has relied on two independent data sets (the NFGS and 
HyperLeda+ samples) and uses a third survey to cali- 
brate completeness (the NYU VAGC low-redshift sam- 
ple). Key results are summarized here. 

• We identify three stellar mass scales of inter- 
est: the "shutdown mass" log M s /Mq ~ 11.2 
above which blue-sequence E/SOs (and the blue 
sequence) generally do not exist; the "thresh- 
old mass" log Mt/M Q ~ 9.7 below which blue- 
sequence E/SOs become extremely abundant; and 
the "bimodality mass" log Mb/M@ ~ 10.5 marking 
the midpoint of this range. Blue-sequence E/SOs 
increase in numbers down to Mb, then level off at 
5-7% of all E/SOs from M b to M t , and finally spike 
up to -20-30% of all E/SOs below M t . Our sam- 
ples do not probe below M* ~ 10 8 Mq, and statis- 
tical analysis of demographics and environments is 
limited to M* > 10 9 M Q . 

• The bimodality mass Mb and the threshold mass 
Mt have been previously linked with strong shifts 
in star formatio n hist ory a nd gas rich n ess b y 
iKauffmann et all (|2003D and IKannappanl (|2004f ). 



respectively (see also iKannappan fc Wei [2008). 
However, our definition of M t corrects an erroneous 
associ ation of M t with Mb made by IKannappanl 
(|2004T ) , caused by assuming agree ment between the 
stellar mass zero point s used bv IKauffmann et al.l 
(|2003ft and iBell et all (|2003l ). The stellar mass 
zero point used in the present wor k agre es in 
normalization with IKauffmann et alj ((2003) and 
IKannappan fc Gawiserl ( 2007h . 



• Morphologies and numbers of blue-sequence E/SOs 
with masses near M s are consistent with merger 
remnants that will fade up to the red sequence 
( "quenching mergers" ) . 

• Blue-sequence E/SOs below Mb have fairly settled 
morphologies and numbers in the population ex- 
ceeding the rate of close pairs, suggesting disk 
(re)building may be common. Most are SOs (in- 
cluding SO/as), and all of the Es have some type of 
disk. 

• Below Mb, we measure abruptly lower density envi- 
ronments for both red- and blue-sequence E/SOs in 
our data sets, often even lower than for many late- 
type galaxies. While these data are incomplete, 
we conclude that the E/SOs in our analysis repre- 
sent a distinct population from dEs, which would 
be found primarily in a lower mass range in cluster 
environments. 

• Blue-sequence E/SOs are intermediate between red- 
sequence E/SOs and spirals in M»-er and M*- 
radius scaling relations, but are more similar to 
red-sequence E/SOs, validating their morphologi- 
cal classification. Given that most sub-Aff, E/SOs 
in our samples reside in modest density environ- 
ments, their compact structure is most naturally 
explained by past merging. We note the emergence 
of a shifted M*-radius relation below — M t . 

• Direct evidence for disk-building in blue-sequence 
E/SOs includes blue outer disk colors and a high fre- 
quency of kincmatically distinct (polar or counter- 
rotating) secondary stellar disks. Statistically, the 
strongest color difference between red- and blue- 
sequence E/SOs is in outer disk color. 

• Based on current gas content and specific star for- 
mation rates, the growth potential of blue-sequence 
E/SOs is comparable to that of late-type galaxies 
below Mb- Moreover, several lines of evidence sug- 
gest that star formation in blue-sequence E/SOs is 
intermittent, involving both minor and significant 
starbursts. We also find evidence that gas and stel- 
lar disk growth may be partly decoupled. 

• Pseudobulges may form in tandem with extended 
disks in many blue-sequence E/SOs. All sub- 
Mb blue-sequence E/SOs in the NFGS have mod- 
est central velocity dispersions, and about half 
show evidence of minor disturbances and centrally 
enhanced star formation, suggesting gas inflow 
events. However, these results are preliminary 
given the small sample size. 
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• Red-sequence E/SOs may also participate in disk 
regrowth. Below Mb, red-sequence E/SOs show 
substantial diversity in radii, central velocity dis- 
persions, gas fractions, and specific star formation 
rates. A subset of these systems could represent 
"off" states in a bursty cycle of disk regrowth. 

• A significant fraction of blue-sequence E/SOs seem 
to occupy a "sweet spot" in parameter space, with 
masses and environments in a regime characterized 
by abundant gas, combined with concentration in- 
dices near the value identified bv lKauffmann et al.l 
(2006) as optimal for peak star formation efficiency. 
This raises the intriguing possibility that the pro- 
cesses that form blue-sequence E/SOs actually en- 
able galaxies to build large spiral disks and escape 
from the dwarf irregular regime. Based on down- 
sizing considerations, we conjecture that some frac- 
tion of z = spirals above Mb may have grown from 
analogous blue-sequence E/SOs more massive than 
those seen today. 



Several areas remain for future work. Spatially re- 
solved analysis of gas reservoirs, star formation, and 
star formation histories is underway to quantify disk and 
bulge growth in both red- and blue-sequence E/SOs. This 
analysis will assess E/SO evolutionary trajectories and 
examine the ratio of fading-merger vs. disk-building sys- 
tems as a function of mass. Evaluating the spiral-galaxy 
progenitor scenario will require comparative analysis of 
high-redshift blue-sequence E/SOs and presumed spiral 
descendents today. A much more extensive analysis of 
the internal dynamics of today's E/SOs with masses be- 
low Mb is also essential to establish the plausibility of 
their building dominant dynamically cold components, 
i.e., pseudobulges and spiral disks. If this picture proves 
correct, then blue-sequence E/SOs may represent the first 
direct evidence for a substantial and well-defined class of 
merger-formed galaxies caught in the act of rebuilding 
disks, as predicted by hierarchical models of galaxy for- 
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